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CHAPTER 



1 



INTRODUCTION 



The radio index of refraction in the atmosphere 



is defined as the ratio of the speed of light to the velocity 
of a radio wave in the atmosphere. There are two methods 
of measuring the index of refraction, the direct method em- 
ploying radio refractometer s and the indirect method utiliz- 
ing standard weather observations. The indirect method re- 
quires the measurement of temperature, pressure, and humidity 
with the subsequent conversion to the index of refraction 
using the formula below: 



(n-1) x 10 6 = 77 - 6 P 



T 



lie + 5.8 X 10 s 
T T 2 



( Bean et al, 1966) 



or a more convenient expression in N units 



N - (n-1) x 10 6 = lh2Z 



T 



lie + 5.8 X 10 5 
T T 2 



where, n = the index of refraction of the atmosphere 
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2 



p = barometric pressure in millibars 
T = temperature in degrees Kelvin 
e = water vapor pressure in millibars 

The operation of radio refractometers, used for 
direct measurement; is based on the fact that the resonant 
frequency of a cavity resonator is inversely proportional 
to the index of refraction of the gas filling the cavity. 
Most ref ractometers utilize two resonant cavities; a sealed 
reference cavity and a ventilated sampling cavity. Since 
the reference cavity is sealed and maintains a constant res- 
onant frequency (within the limits of small temperature 
and pressure changes) ; the change in resonant frequency of 
the sampling cavity due to the change in the index of re- 
fraction of the air passing through it can be detected. 

There are a number of ways of detecting this resonant fre- 
quency change and converting it to the index of refraction; 
for example; the Crain type r ef ractome t er measures the dif- 
ference in frequencies and; 



Af 

f 



-An ; AN = An X 10 6 



The Birnbaum type ref ractome ter measures the phase differ- 



ence between the two cavities. 
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Considering only the accuracy of the measurement 
and the time required to evaluate the data, the direct 
method is desirable since the value of N is determined 
directly and only one sensor is required instead of three. 
However, the disadvantages of high cost, complexity, and 
maintenance required outweigh the advantages and the majority 
of refractive index measurements to date have been made in- 
directly. 

The ideal situation, especially in the field of 
radio propagation studies where precise measurements of the 
index of refraction are required, would be to have a simple 
inexpensive method of measuring the index of refraction di- 
rectly. One possible method of accomplishing this is to use 
the backscattering characteristics of a ventilated resonant 
cavity as the cavity falls through the atmosphere. As the 
cavity falls through the atmosphere, it can be tracked by 
radar and the change in backscattering due to a change in 
the resonant frequency can be detected and related to the 
index of refraction of the gas passing through the cavity. 

In recent years a considerable amount of research 
has been directed to the reduction of backscattering or 
radar cross-section by impedance loading particularly in 
the resonant region. Chen and Liepa (1964) studied the 



4 



scattering cross-section of a thin cent er -loaded cylinder 
and demonstrated the reduction of the radar cross-section 
due to loading, Liepa and Senior (1966) studied the scat- 
tering behavior of a metallic sphere with circumferential 
slot loading. Plonus (1967) investigated the scattering 
behavior of a spherical shell with a finite slot for loading. 
Chang and Senior (1969) studied the scattering behavior of 
a thin spherical shell with a circular aperture and demon- 
strated that there is a large reduction of the radar cross- 
section in the vicinity of the cavity resonant frequency. 

This thesis is concerned with determining if the 
resonant cavity approach to the direct measurement of the 
radio index of refraction in the atmosphere is feasible. 

A thin spherical copper shell with circular apertures has 
been selected for this study because of the symmetry pro- 
perties of the sphere and the circle. 



C H A P T E R 



I I 



THEORY 



2.1 Solutions to Wave Equation in Spherical 
Coordinates with Axial Symmetry 

With axial symmetry Ramo, Whinnery and Van Duzer 



(1965) approaches the problem by separating the solutions 
into waves with components E r , Eg, and H<j> and waves with 
components H r , Hq, and E <j> . Using the spherical surface 
(r = constant) as the transverse surface the wave with 
components E r , Eq, and H is called the TM type and the 
one with components H r , Hq, and E 4 > is called the TE type. 



E^ = 0 and S/S<t> = 0, Maxwell's equations in spherical 
coordinates become: 



Considering only the TM type with H = 0, H = 0, 

T (y 



VX H = j(B e E 



r s in 0 



a 

r 



S ( H sin 0 ) 



S0 




V X E 



jcupH = a ^ 



S(rE ) Se 
0 



r 



( 1 - 2 ) 



Sr S © 



5 



6 



then since H — a^H^ and E = a r E r + agEg 



iH. . „ 

r Sr r S© 



(1-3) 



r s in 6 



S ( sin 9) 

S© 



= jcd e E 3 



(1-4) 



_! ^ rH ») 



r Sr = £E 0 



(1-5) 



solving for E^ and Eg in equations (l-4) and (1-5) we get 



E = 



E , = 



S(H $ sin 0) 



jco e r sin © S© 

1 MrH $ ) 



0 j to £ r Sr 



( 1 - 6 ) 



(1-7) 



substituting equations (1-6) and (l-7) into equation (1-3) 
and carrying out the required operations, we get: 



* (*H,) 

Sr ^ 



+ ! 



1 S ( H (j, sin ©) 



sin © 



S© 



S© 



2 

+ co p. 



e rH* = 0 



( 1 - 8 ) 



9 2. 

Letting k c = co p.e (Ramo et al., 1965) solved the second 
order differential equation (l-8) by the product solution 



method and found: 
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H* = 



f V(V3) (kr) + 7=V(1/Z) (kr) 



p (cos e) 
n 



(1-9) 



where, Jn+(l/2) = a Bessel function of the first kind 
N n +(i/2) = a Bessel function of the second kind 
P^ = a Legendre polynomial 

If the problem contains the origin B n must equal zero be- 
cause W n + (]_y'2) ■'• s infinite at the origin. If the problem 

extends to infinity the Bessel functions must be combined 

( 2 ) 

into the second Hankel function, H n +(]_/2) = ^n+(l/2) “ 

^n+(l/2) re P resen t the wave traveling away. Ramo et al, 
(1965) uses the notation Z n+ (i/2) to represent the particu- 
lar combination of J n +(i/2) and ^n+(l/2) required for the 
problem under consideration and then: 



H * = — P 1 (cos 9 ) Z , . . (kr) 

* n' n+(l/2) 



( 1 - 10 ) 



Substituting equation (1-10) into equations (l-6) and (1-7 ) 
and differentiating, we get: 



A P^( cos 9) 
t? _ n n 

E e - ..3/2 



jcu e r 



nZ n + (l/2) (kr) 



krZ n-(l/2) (kr) 



( 1 - 11 ) 



8 



- A nZ , , . (kr) 

E = n n tU J U 



jco e r 



3/2 



s in 0 



cos 9 pjj(cos 0) - P^ + ]_(cos 0) 



( 1 - 12 ) 



The same approach can be used to determine the 
field equations for the TE mode where b/bb = 0, E r = 0, 

E© = 0, and H<t> = 0 or by the principal of duality replac- 
ing E r by H r , E© by Hq, H<j, by - E^, e by n, and using a 
new constant B n in equations (1-10), (l-ll), and (1-12). 

Applying the principal of duality the following equations 
are obtained for the fields of the TE mode: 



B 

E <t> = -p_P„(cos ©)Z n + (i/2) (kr ) 

nZ n+(l/2) ( kr ) " krZ n-(l/2) ( kr ) 



-B B (cos yj 



u _ n n 
H© 

.i cDur 



3/2 



(1-13) 



(1-14) 



B nZ ... (kr) 
H v = _S n+ ( 1 / 2 ) 






3/2 



sin 0 



cos 9 P^(cos 0) - Pn+q(cos 9) 



(1-15) 



2.2 Spherical Cavities 

The TE-^q^ and modes of a spherical cavity 

will be analyzed using equations (1-10) through (1-15). 
Assume the cavity to be closed with inner radius a and 



that the shell is perfectly conducting. The Bessel func- 
tion of the second kind must be eliminated from our equa- 
tions since the origin is included in the area of interest; 
therefore, the field equations for the TEqoi mode in the 
cavity are: 




( 2 - 1 ) 



-b p 1 (cos e) r 

11 




kr J i/2 ( kr ) 



( 2 - 2 ) 





cos 0 P^(cos 0) - P^(cos 0) (2-3) 



where P^(cos 0) = sin 0 

Pi (cos 0) = 3 sin 0 cos 0 

U 





sin kr 



Making the substitutions indicated above in 



equations (2-l) through (2-3) we get the following equations 
for the mode in the cavity: 
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E* = -L 



[~2 


"s i n k r "1 


- 


- cos kr 


b tk 


kr J 



(2-4) 



H e = 



■B^ sin 9/2 



jcoiar 




- kr sin kr - cos kr 



LVkr 



(2-5) 



H r = 



jo»iir 





sin kr 


- 


- cos kr 


Jnk 


. kr 



( 2 - 6 ) 



Since we made the assumption that the shell of the cavity 
was perfectly conducting, the tangential electric field at 
r = a must go to zero, therefore; 



E 



4 > 



B s in 0 


fz 


"sin ka 


n ^ .1 , 


1 


— n n ° V ^ 


a J 


rtk 


ka 



0 = - S J- - A . . ^ a - cos ka 
ka 

0 = tan ka - ka 
tan ka = ka 

Tan ka = ka is the condition for resonance for 
the TEj_oi mode in a perfectly conducting closed cavity. 

The solution of this transcendental equation can be solved 
graphically and the solution is ka~4.5. 

For the TM^q-^ mode we have: 



11 



H* = 



A^ sin 9 



2 


sin kr 


1 


- cos kr 


Ink 


kr 



(2-7) 



E e = 



A^ sin 9 
2 

jcn e r 




— kr sin kr - cos kr 



LKr 



( 2 - 8 ) 



E r 



2A cos 9 I 2 

1 

. _ 2 
jo) e r 



sin kr 



kr 



- cos kr 



(2-9) 



As for the TE^q^ mode the tangential electric 
field must equal zero at the shell for the resonant condi- 
tion for the TM^oi mode, therefore: 



0 = sin k * JL - ka - 1 
cos ka ka 



tan ka 



1 




ka 

1 - ( ka) ^ 



The graphical solution to this equation is 

k a* 2 . 7 5. 

The figure of merit or ,! Q," of the cavity is de- 
fined as the ratio of the energy stored at resonance to 
the energy loss. Ramo et al. (1965) defines the Q for the 
TMioi mode as approximately equal to the ratio of the in- 
trinsic impedance of the gas filling the cavity to the 
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surface resistivity of the metal in the shell. Then for a 
copper cavity filled with free space the Q, would be; 



Q, Ra 



R 

s 



377 



- 7 

2.61 x 10 




for f = 10 GHz 



3 



Q 

Q 

Q 



^ 14.4 x 10 

is sometimes defined as 
f 

f^S 

Af 



where f = resonant frequency of the cavity 
Af = change in frequency between half 
power points 

Then the Af for this mode is approximately equal 
to 695 KHz. 



2.3 Scattering by a Spherical Shell 

To study the scattering behavior of a spherical 
shell, it will be necessary to expand the solution of the 
wave equation in section I to include variations with 0. 
This expression or complete solution of the wave equation 
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in spherical coordinates is given in terms of spherical 
vector wave functions by Stratton (1941). There is no 
loss in generality to assume that we have a plane wave 
incident on the sphere in the negative z direction with 
the electric field vector in the positive x direction, in 
a spherical coordinate system ( r, 9, <i>) related to the 

Cartesian coordinate system by x = r sin 9 cos <t>, y = 
r sin 6 sin ( t ) , and z = r cos 9. Making these assumptions 
the incident field equations can then be written as: 



E 

H 




jkz 

e 




jkz 

e 



j W "b 

Normalizing the E vector to unity and suppressing the e 
we have: 



E = a e^ kZ (3-1) 

x 

- - jkz 

H = -a Y e (3-2) 

y o 



where k = propagation constant 

Z q = intrinsic impedance of free space 

Y = intrinsic admittance of free space 
o 
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Employing the spherical vector wave functions 
mentioned above, the incident field is expressed in spheri- 
cal coordinates as: 



i V — ^ n 

E = X j 2n + ] 

Z, "(n + 



1 ) 



n=l 



1 v — < n 

H = jY \ j 2n + 

° / n ( n 



+ 1 ) 



(M^ - j ) (3-3) 

v Oln eln' v ' 



( 1 ^ 1 ^ ) ( 3- 4) 

Oln eln 



n=l 



where and ^ are the spherical vector 

wave functions as defined by Stratton (1941) as given below 



m! 11 

gmn 


1 + 
II 


m f n( kr) 
kr 


< 


cos 0 ) 
sin (9 


s in 

m<j> 

cos 


a e 




- 


♦n(kr) 

kr 


( 


cos e )] 

~5e 


cos met) 
sin 


a <J) 


(1) 

N e 

§mn 


= n( n + 1 ) 


»n(.kr) 
(kr) * 


m 

. P n ( 


cos e) 


COS 

m<t> a 

. r* 

sin 




+ 


^n(kr) 

kr 


d [ P n (cos $ )J 

Se 


1 cos 

sin 


m 4> a 

G 




+ 


B. ..*.n.(£r 
kr 


1 p n (cos 9) 
sin Q 


sin 

m<t» a ^ 

cos 
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where 

\|< n ( kr ) = kr J n (kr), J n (kr) is the spherical 
Bessel function of order n and the 
prime indicates differentiation with 
respect to the entire argument, 
m 

P (cos d) = A Legendre polynomial 
n 

Substituting equations (3-5) and (3-6) into 
equations (3-3) and (3-4) the incident fields are given 
as follows: 



E 



i 



I 



n 



2n + l 
n(n + 1) 



Ukr) 

kr 



■ n (cos g) cos $ 
sin e 



- ^ n ( kr) 
kr 



c)D£ (cos 9)1 

*0 



s in 0 a , 
4> 



-j n( n + 1) tn(kr) (cos 9) cos <j> a 

(kr)* 




LMl 



kr 



3 P*n (cos 9 )] 



cos 






a 

9 



+0 tyh (kr) 
(kr) 



■ n 



(cos e ) 
sin e 



sin <J> a , 
<t> 



(3-7) 
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i , n 

H = JY 



J 2n + 1 
n(n + 1) 



n=l 



n(n + 1) ^n( kr) P n (cos 0) sin $ 

kr 



+ ^n( kr ) 


d[P n (cos 8)] sin 4> 


kr 


^0 


^n(kr) 


P n (cos 9) cos <t>a A 


kr 


sin Q 


+j tn(kr) 


P n ( cos 0) s in tj) a 


kr 


sin 0 


+d ^n(kr) 


d L P rj ( c o s aJ cos 


kr 


S0 



(3-8) 



Then breaking the incident fields into their 
respective components, we have: 



E = -d 



00 n 

cos 4> X ^ ,j n ( 2n + 1) (kr ) p n (cos 0) (3-9) 



(kr) 



I 



n=l 



kr 



E a = 



COS 0 

kr 



00 

2n + 1 
/ . n(n+l) 



^n(kr) Pn( cos 0 ) -j ^n(kr) ( cos 0)] 



sin 6 



n=l 



(3-10) 



E = sin 4> 



4> 



kr 



^ n 

\j Sn + 1 : 
/ m n ( n +1) 



3 fn(kr) Pnfcos 0) -ty n (kr) d[ P n ( c o s g)] 



sin e 



$9 

(3-11) 



n=l 
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H 



l _ 



00 

Z' 

n=l 



jY sin 0 j ( 2n + 1) t n (kr) P„(cos 0) 

° (kO 2 



H = jY sin $ 
0 o 

kr 



2n ♦ 1 

/ . n ( n + 1) 
n=l 



'l'n(kr) 3 [ p n (cos 0)] + 



J ’•'n(kr) P n (cos 0) 
sin 0 



H = 



jY cos <J> 
o 

kr 



2n + 1 0 ^n(kr) 5 [ P jj( cos 0 )] 

n(n + 1) ci0 



n=l 



^n(kr) P j (cos 0) 
sin G 



( 



The scattered field can be represented in the 
manner as the incident field with spherical vector wave 
functions as follows: 



(4) 



(4) 



E = > (A m M . - 0 B„ N , ) 

n oln n eln 



( 



(4) 



(4) 



H ^oln ^ p n ^eln 

n=l 



3-12) 



3-13 ) 



3-14) 

same 



3-15) 



0 



(3-16) 
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Where ^ and ^ are the same as ^ and except 

that to satisfy the requirement that the scattered field 

represent an outgoing field at infinity, i n (kr) is replaced 

by p (kr) = kr \kr) and H ^ c ^ is the spherical Hankel 
n a n 

function of the second kind and A and B are unknown 



amplitude coefficients. 

The boundary conditions that the tangential elec- 
tric fields at the surface of the sphere must equal zero 
provide a means for solving for the unknown amplitude coef- 
ficients A and B as follows: 



n n 



E t = E 1 + E S 

at r = a 

E 1 + E s = 0 

e e 

E 1 + E s = 0 

<J> <t> 



then 



t 

E 

0 



+ 





CO 

\ ' n 




0 = cos 4> 


> d 2n + 


1 


ka 


n( n h 


- 1) 




n=l 






! 

0 +n 


(ka) 


00 

cos (a 


P n ( k a ) Pn 


(cos 


- 21 " 




sin i 



n=l 



Mka) P n (cos 0) - (3-17) 

sin 0 

d Pn (cos 0 )] 

^0 

el -d B n Pn (ka) 3 [p P ( cos 6\ 1 

^0 



E t = 0 



sin (i) 
ka 




n 



j 2n + 1 
n( n + 1 ) 



j (ka) ?n (cos 9 ) 
sin 0 
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(3-18) 



(ka) 5 [ p n( cos 9 )3 



00 

sin <j> % p B p' (ka) P n (cos 0) 

ka n sin Q 

n=l L 



A p ( ka ) d [ p n ( cos 0)] 



n n 



Solving equations (3-17) and (3-18) for and 



B we get: 
n 



n 



A n = -j 2n + 1 



\ n ( ka ) 



n(n + 1) f> n (ka) 



.n 



B n = -J 2n + 1 



(ka) 

n(n + 1) pj* (ka) 



(3-19) 



(3- 20) 



The scattered transverse electric field components 
for large kr can be approximated by replacing p^(kr) and 
(kr) by the leading terms of their asymptotic forms 
as shown by Senior and Goodrich (1964) and the results are: 

1 



„s . - jkr 

E = j cos <t> e 

kr / . n( n + 1 ) 



r\^ .n „ 

\ j 2n + 1 

n( n + 1 

n=l 



^n (ka) 5 [ p n ( cos 9 )] 



L Pn ( ka ) 



Se 



(3-21) 



^n(ka) P n ( c o s 9 ) 



P n (ka) sin 0 _ 
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= - j sin (t> e 
kr 




n=l 



j n ?n + 1 
n(n + 1) 



in(ka-) P^(cos 0 ) 

p ( ka ) s in 0 
n 



t n ( ka ) d [P^( cos g)] 

P n (ka) Se 



(3-22) 



E, 



= 0 



A more concise way of writing the equations of the scattered 
fields is: 

E S = j cos d> e J * kr S S ( 0 ) (3-23) 

9 ” 1 
kr 

o -jkr ~ 

E = -j sin d) e 3 (0) (3-24) 

* kr 2 

S s 

where: S^(e) and S^^) are the summations in 

equations (3-21) and (3-22). 



2.4 Scattering by a Spherical Shell With 
a Circular Aperture 

This problem is separate from the one in the 
preceeding section, but it is related to it. As we have 
shown, the incident and scattered fields can be represented 
in spherical vector wave functions as follows: 



21 



E = 



2 

n=l 



n ,(•!■) . ( 1 ) , 

0 to + 1 («oin - ^ B el„> 

n ( n + 1 ) 



H 



00 

2- 

n=l 



n . ( 1 ) ( 1 ) 

= ,j Y > j 2n + 1 ( N _ . - .1 M . 

o / — r Oln e in 

n( n + 1 ) 



(4-1) 



(4-2) 



E = 



00 

2 

n=l 



(4) (4) 

( -^n ^“Oln “ j ®eln^ 



H‘ 



00 

2 

n=l 



(4) . (41, 

- y (-^n ^Oln “ 0 ^eln) 



(4-3) 



(4-4) 



And applying the same technique the fields inside the shell 



can he represented as: 

00 



2 r U) • U \ 

n ^Oln " ^n ^eln] 



n=l 



H 



Vr (1) • (1) 1 

in - Y o / 1? n N 03 n ^ D n M elnj 

AMffP 



(4-5) 



(4-6) 



n=l 



The far field scattering coefficients are now: Chang 

and Senior (1969) 



( 9 ) = j J"A n Pn 

n=l L 



(cos e ) - B n d[P n (cos 9 )] 

^0 



sin 9 



5 z («) - 



00 

= / j n A n a Oji (cos e )] - B n Pn( COS 6 ) 

sin 9 

n=l L J 



(4-7) 



(4-8) 
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And: 

E® = -j sin <t> exp( - ,jkr ) (e) 

kr 

E® = j cos 4> exp( - ,-jkr ) S s (e) 

G -L 

kr 

And from Senior and Goodrich (1964): 
cr(0, '*’) = (T e (e, < t > ) + cr <(>( 0 , ) 



where 



o- 0 ( 9 , I* ) 



2 2 
cos ^ 4> 

n 



sj (0) 



♦) 




. 2 
s in 





(4-9) 



(4-10) 



(4-11) 



( 4-12) 



(4-13) 



Equations (4-1) through (4-8) show that once 

the amplitude coefficients A , B , C , and D are deter- 

n n* n n 

mined the fields everywhere are known. Therefore, the 
problem is a boundary value problem where the boundary 
conditions are the tangential electric fields are zero 
on the surface of the shell and the fields through the 
aperture must be continuous. When these boundary condi- 
tions are applied to the equations, it is possible to 



solve for the coefficients. 
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Chang and Senior (1969) discuss a number of ap- 
proaches to this problem for example, using integral equa- 
tions in terms of the unknown aperture field, and the 
variational technique which relates the backscattering 
cross section to a scalar function U. U is in terms of 
the incident electric field and the total surface current 
density of the shell. Both of these approaches require 
the choice of a set of expansion functions; for the integral 
equation approach a set of expansion functions is chosen to 
represent the aperture field and for the variational method 
a set of expansion functions is chosen to represent the 
surface current density. The accuracy of the solution for 
both cases is determined by the proper choice of these 
expansion functions. Due to the presence of the cavity 
in this problem the proper choice would be extremely dif- 
ficult . 

Due to the difficulties mentioned above, Chang 
and Senior (1969) applied the method of the least square 
error and then modified this method to separate out the 
discontinuties at the edges of the aperture. This ap- 
proach will now be summarized: (Refer to figure 2- 4a 

for the geometry for this solution). 
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Figure 2-4a. Coordinate System 
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Applying the boundary conditions to equations 
(4-1) through (4-6) and solving for the coefficients A 



and B in terms of C and D we get: 
n n n 



A p ( ka ) = 
n n 



; n 2n + 1 



C - j 
n n( n + 1 ) 



t n ( ka ) 



(4-14) 



B n p n (ka) = 



D - 2n + - 1 ■ 

n n( n + 1) 



( ka ) 

n 



( 4-15) 



Making these subsitutions, the equations for the scattering 



coefficients A and B are: 
n n 




A p (ka) n (cos 0 ) b p ( ka ) 5[Pn( cos 0 jl 

n n . 0 n n v ' L - - ■ - - s 

sin 9 oq 

K'L P (ka) 1 £Ll££^£)L 3 b c ' (k a) P n 
/ J n n W n n sm 9 



( 4-16) 
= G 1 (0) 

( 4-17) 
" G 2 (0) 



n=l 



for 9 < 9 < n 

o ~ 

And; 



00 

l 

a 



rA n ^[ P n (° os 0 )J _ L i B n p n ( cos 

^ n (ka) §0 \]; n '(ka) sin 9 



= 0 (4-18) 



n r n (cos 0) B n S [ P i ( cos ^ )] 



sin 0 



^n ( ka ) 



^0 



= 0 (4-19) 



for 0 < 0 < 0 
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where 



G 1 (0) = - 



00 

V i n 2n + : 

Z ■>(" * 



1) 



n=l 



* (ka) F n (cos a) 



n 



sin e 



-j * (ka) 



dfn (cos 0 )J 
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(4-20) 



G„(0) = -> j n 2nj_l 



n( n + 1 ) 

Pn ( cos e) 

0 t n ( a ) gj_ n Q 



ty n ( ka ) ^ & n (^ 0S ,_g_)] 



(4-21) 



Applying the method of least square error, Chang 
and Senior (1969) obtained a set of 2M X 2M matrix equations 
which when solved give an approximation for the 2M scat- 
tering coefficients A and B„ but due to the discontinuities 

n n 

at the edge of the aperture these equations converge very 
slowly. To improve the solution Chang and Senior (1969) 
examined equations (4-16) through (4-19) and separated out 
the discontinuities. Then applying the least square error 
method to the new equations they arrived at a linear system 
of 2(M + 1) equations with unknowns defined by: 



l n = a n + K l * n (ka) 7 n 



B„ = + K 



n - n 1V 2 -r—i 



P n '(Ka) n 
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where 





2 



d 



n ( 2n + 5)( 2n + 3 ) dn+2 



c 



n ( 2n - 1)( 2n + 3) 



1 



( 2n - 3 ) ( 2n + l) 



d 



n 



exp [j( n + l/2)e o ] 



which converge more rapidly giving a more accurate approxi- 
mation for the scattering amplitude coefficients for a 
given M than the original least square method. 

The results of this "modified least square method" 
produce the following for far field scattering coefficients: 
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(4-23) 



X 



2n ,i,i 



(ka) n 



X 2M + 1 = K i 
X 2M+ 2 K 2 



The complete development of this approach and 
the computer programs to carry out the calculations are 
presented in Chang and Senior (1969). The theoretical 
response for an aperture size of 0 Q = 30° with 0=0, 

<t> = 0, and 0.8 < ka < 4.8 with M = 10 for ka < 1.1, M = 20 
for 1.1 < ka < 4.8 and M = 30 in the region of expected 
peaks as calculated by Chang and Senior (1969) is shown 
in figure 2-4b. 



10.0 
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Figure 2-4b. Computed Backscattering Cross Sections of a Spherical Shell with 
0 Q = 30° (Using the Modified Least Square Method). 
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2.5 Theoretical Case of A Remote Cavity 

The theory in the preceeding sections will now he 
applied to a remote cavity. The assumption is made that 
the Q of the cavity with an aperture is the same as the 
closed cavity. It is known that the aperture in the cavity 
will lower the Q, hut the calculations employing the per- 
turbation theory as in Slater (1950) are extremely compli- 
cated even for a small perturbation and are beyond the 
scope of this thesis. 

The fact that the back scattering is reduced in 

the vicinity of resonance will now be employed to determine 

the detectable change in frequency and subsequently the 

change in N units. Assume that the resonance frequency of 

the cavity is 10 GHz (a = 2.15 cm). 

For a change of 100 N units, Af = f - jr.Q. = 999901 

0 1.0001 

Or Af for a change of one N unit is approximately 
10 KHz. 

Assuming that a change of one half power can be detected 

accurately leads to the assumption that a change of one 

half of the frequency between half power points can be 

detected. Using the theoretical maximum Q for the ™j_oi 

mode from section II of 1.44 x 10^ we get = 347.5 KHz 

2 

and represents a change of approximately 35 N units. 



CHAPTER 



III 



EXPERIMENTAL RESULTS 



To confirm the theory in the preceding sections, 
an experimental study was conducted. 

3.1 The Model and The Experimental Set-up 

A hollow spherical copper shell with a 2 inch out- 
side diameter and a single circular aperture of 9 0 = 15°, 
Figure 2-4a, was employed as the resonant cavity model to 
investigate backscatter ing near resonance from such a cavity. 
Due to the shape of the field lines coupled to the inside of 
the resonant cavity for various modes, the TE, mode was 
selected because the circular aperture causes less distor- 
tion of the field lines for this mode. The theoretical free 
space (vacuum) resonant frequency calculated using the 
outside diameter and neglecting the aperture was 8.45 GHz. 
(Found later, experimentally, to be 8.55 GHz.) 

An anechoic chamber 6 feet wide, 8 feet high and 9 
feet long shown in Figure 3-la was constructed to eliminate 
extraneous reflections and outside interference. The resonant 
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Figure 3-la. Photograph of Experimental Setup 
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spherical cavity model to be examined was placed on a 
styrofoam pedestal inside the chamber three feet from the 
illuminating antenna. 

The equipment arrangement for measuring the 
backscattering is indicated in Figure 3-lb. The HP 624B 
signal generator was used as the transmitter and the TS 
147D signal generator was employed as the local oscillator. 
By frequency modulating the TS 147D with the sweep genera- 
tor, the requirement for continuous frequency tracking of 
the transmitter was eliminated. Magic tee number two 
served as both the mixer and detector. The 1216A ampli- 
fier was used as an IF amplifier and the output read di- 
rectly off the oscilloscope face. The oscilloscope response 
represented the band-pass characteristics of the IF ampli- 
fier, but the peak response was directly related to the 
amplitude of the electric field reflected from the spher- 
ical shell. 

The procedure for taking a measurement was as 
follows: With the transmitter tuned to the desired fre- 

quency and the spherical cavity model removed from the 
anechoic chamber, arms A and B of the first magic tee were 
balanced by tuning the three tuning screws in arm A and 
the variable attenuator in arm B until no signal was 



34 





S- 




o 


CL 


+-> 


CD 


rtJ 


CD 


S- 




CD 


LT) 


C 




CD 











o' 


O 




X 




CD 




i— 


CL 




•r^ 


O 




O 


O 


>> 


CO 


CO 




o 








i 

CO 

CD 

s- 

3 

CT> 



Block Diagram of Experimental Equipment 



35 



observed on the oscilloscope. This meant that none of 
the transmitter signal from arm C was entering arm D and 
that the transmitter signal was divided equally between 
arms A and B. The spherical cavity model was then placed 
on its pedestal and irridiated. One half of the signal re- 
flected from the model entered arm D and was mixed and 
detected in magic tee number two. The detected signal was 
then amplified and recorded. The model was rotated from 
its initial position of 0 = 0° to 0 = 90° in 10° steps and 
the measurements repeated for each step. The final step 
was to rotate the model to 0 = 180° (circular aperture on 
the back) and make a measurement, 

3.2 Measurements 

The initial measurements were made to determine 
the resonant frequency of the cavity. The copper shell 
has a finite thickness which causes the cavity to be 
smaller than the outside measurement would indicate and 
increases the resonant frequency while the aperture has a 
tendency to lower the resonant frequency. The theory shows 
that the signal backscattered from the aperture (0 = 0°) is 
less than from a solid shell (© = 180°) only in the 
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vicinity of resonance. For the above reasons, measurements 
of the backscatter from the aperture and the solid shell 
were taken starting from a frequency below the theoretical 
resonant frequency and recorded as frequency was increased. 
The experimental resonant frequency was found to be approxi- 
mately 8.55GHz. 

After the resonant frequency was determined, 
measurements of backscatter ing were made in the vicinity 
of this frequency. The ratio of the backscattered signal 
from the shell with 9=0° to the backscattered signal from 
the shell with 0 = 180° was recorded and plotted. The 
first attempts to obtain repeatable data failed and was 
believed to be due to the inaccuracy of the frequency 
measuring device on the transmitting equipment and a 
considerable amount of copper oxide present on the copper 
shell surface. The shell was cleaned and the measurements 
repeated. It remained extremely difficult to obtain re- 
peatable data, but an attempt was made and the results of 
three sets of data were averaged to obtain the curve shown 
in Figure 3-2a. The curve in Figure 3-2a was obtained by 
computing the ratio of the signal backscattered from the 
shell with 9=0° to the signal backscattered from the 
shell with 9 = 180° and the results plotted versus frequency. 
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Measured Backscattered Signal of a Spherical Shell 
with e 0 =15° and 0=0° 



38 



This curve provides a means of determining the approximate 
bandwidth, the resonant frequency, and the detectable 
change in N units. 

As the backseat ter ing measurements above were 
being made, the variations in the backscatter ing as a func- 
tion of 9 were also recorded in 10° steps from 0=0° to 
0 = 90°. The purpose of this measurement was to determine 
the effects of the aperture as a function of illumination 
angle much as what would happen in the use of a free fall- 
ing remote cavity. The results of these measurements are 
shown in Figures 3-2b through 3-2d. These curves were ob- 
tained by taking the ratio of the signal backscatter ed from 
the shell at an angle 0 to the signal backscatter ed from 
the shell at 0 = 180° with the frequency held constant and 
the results plotted as a function of 0. Figures 3-2b 
through 3-2d show the variations in the backscatter ing as 
a function of 0 for three different frequencies, Figure 
3-2b was plotted for a frequency slightly below resonance 
while the data for Figure 3-2d was taken at a frequency 
slightly above resonance and the data for Figure 3-2c was 
taken at the approximate resonant frequency of the 



cavity. 



39 




co 

-o 



Measured Backscattered Signal of a Spherical Shell with 
6o - 15° as a function of 6 for f = 8.549 GHz 



Shell at 0) 

Shell at 0 = 180°) 



40 




CD 

"O 



Measured Backscattered Signal of a Spherical Shell With 
6 0 = 15° as a function of 6 for f = 8.550 GHz 
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3.3 Discussion of Experimental Results 

The objectives of the experiment were to: 

(1) Determine the resonant frequency of the cavity 
with a circular aperture. 

(2) Determine the band width of the cavity with a 
single aperture. 

(3) Determine the change in N units, of the gas fil- 
ling the cavity, that could be reliably detected. 

(4) Determine the effects of angle of incidence on 
the copper sphere relative to the location of 
circular aperture. 

The resonant frequency was found to be approxi- 
mately 8.55 GHz and the backscattering was investigated in 
the vicinity of this resonant frequency. Figure 3-2a re- 
presents the results of this investigation. From the curve 
in Figure 3-2a, the approximate "Q" of the cavity with a 
single aperture was found to be 1600 with a corresponding 
Af of approximately 5.4 MHz. This Af would allow the de- 
tection of a minimum change of approximately 320 N units. 
Since the maximum change in N units in the atmosphere from 
a vacuum to sea level is approximately 340, this system 
would reveal essentially no useful information about the 
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index of refraction of the atmosphere as a function of 
height. To detect a change of one N unit, the "Q" of the 
cavity would have to be 505,000 or have a Af of approxi- 
mately 17KHz at a resonant frequency of 8.55 GHz. 

The most difficult part of the experiment was 
that of obtaining repeatable data. This was due in a 
large part to the inaccuracy of the frequency measuring 
device on the HP 624B. The accuracy of the frequency 
measuring device was ± 0.03 percent or approximately 
_+ 2.565 MHz and this error could be as much as one half Af, 
the bandwidth. This represents a significant error, how- 
ever, the frequency measurements were accurate enough to 
definitely show that the bandwidth of the resonant cavity 
was much greater than the bandwidth needed to be able to 
detect a change of one N unit in the dielectric filling the 
cavity. Another error was introduced in the measurements 
by a noise level in the measuring system of 0.4V. This 
created a problem only when the desired signal level 
dropped to this value and really had very little effect 
on the measurements. 

To satisfy the fourth objective of the experiment, 
the measurements for Figures 3-2b through 3-2d were taken. 
Investigating Figure 3-2c, we see that there is little 
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effect caused by a change in the illumination angle 0 of 
± 20 degrees relative to the location of the circular aper- 
ture but a considerable effect for changes in 0 of more 
than 20 degrees. Figure 3-2b shows that the changes of 
the backscattering with changes of 0 are essentially the 
same when the frequency of the probing signal is slightly 
below resonance as when it is at resonance. However, when 
the frequency is slightly above resonance the changes are 
very different from those at resonance. 

Measurement errors that could be introduced by 
inaccurate positioning of the model are considerable, for 
example, if the spherical copper shell had been rotated 
to an angle of 0 ± 25°, Figure 3-2d, the resonant fre- 
quency would appear to be 2 MHz higher than the actual 
resonant frequency which would represent a change of 234 N 
units in the dielectric filling the cavity. On the other 
hand, if the shell had rotated more than ± 60° there would 
be no change in the backscattering for small changes in 



frequency . 



CHAPTER IV 
SUMMARY 



The feasibility of remote measuring the index of 
refraction of the atmosphere using the resonance charac- 
teristics of a hollow spherical cavity with apertures was 
investigated both theoretically and experimentally. The 
theoretical study assumed a closed cavity and produced a Q, 
with a bandwidth narrow enough to enable the detection of 
a change in the index of refraction of the gas filling the 
cavity of approximately 35 N units. The experimental 
study produced a minimum detectable change of approximately 
320 N units. 

To produce a feasible system for measuring the 
index of refraction of the atmosphere a minimum detectable 
change of no more than 1 N unit would be required. As 
shown by both the experimental and theoretical studies, 
the sensitivity of the proposed system using a hollow 
spherical resonant cavity with one aperture is far from 
that required. In actual practice a number of apertures 
would be required for adequate flushing of the cavity and 
to make the measurements less dependent on the angle of 
illumination, 0. 
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With more accurate transmitting equipment; the 
measurements could be repeated and plotted more accurately, 
but since the equipment used in this experiment show con- 
clusively that the bandwidth of the cavity is far too 
large to enable the measurement of a minimum value of N 
with only one aperture, it is not recommended that this 
method of measuring the index of refraction be further 
explored . 
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